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ABSTRACT 

Accretion onto black holes in wind-fed binaries and in collapsars forms small rotating 
discs with peculiar properties. Such "mini-discs" accrete on the free-fall time without 
help of viscosity and nevertheless can have a high radiative efficiency. The inviscid 
mini-disc model was previously constructed for a non-rotating black hole. We extend 
the model to the case of a spinning black hole, calculate the structure and radiative 
efficiency of the disc and find their dependence on the black hole spin. If the angular 
momenta of the disc and the black hole are anti-aligned, a hydrodynamic analog of 
Penrose process takes place. 

Key words: accretion, accretion discs — black hole physics. 



> 
o 
in 

On 

O 

On 
O 



1 INTRODUCTION 

The mini-disc model was motivated by the estimate for an- 
gular momentum of accretion flows, /, in wind- fed X-ray 
binaries (Illarionov & Sunyaev 1975; Shapiro & Lightman 
1976; Illarionov & Beloborodov 2001). These systems hap- 
pen to have / ^ r^c, where Vg — 2GM/c^ is the gravitational 
radius of the black hole. It is marginally sufficient to form 
a centrifugally supported disc. Then a small disc can form, 
which is not supported centrifugally and instead accretes on 
the free-fall time (Beloborodov & Illarionov 2001, hereafter 
BIOl). The mini-disc accretes so fast (super-sonically) that 
the effects of viscosity can be neglected. A similar disc may 
form inside collapsing stars (Lee & Ramirez- Ruiz 2006; Be- 
loborodov 2008). 

The mini-disc can be thought of as a caustic in the 
equatorial plane of a rotating accretion flow. It absorbs the 
feeding infall, and this interaction releases energy, making 
the accretion radiatively efficient. With increasing angular 
momentum, the size of the disc grows up to lAvgC^ and at 
this point the centrifugal barrier stops accretion, so that it 
can proceed only on a viscous timescale. Thus, the mini- 
disc model fills the gap between two classical regimes of 
accretion — spherical (/ < VgC^ Bondi 1952) and standard 
accretion disc (/ ^ VgC, Shakura & Sunyaev 1973) — and is 
qualitatively different from both. 

The calculations of BIOl were limited to the case of a 
Schwarzschild black hole. In the present paper we study the 
mini-disc around Kerr black holes. The model is constructed 
under the following assumptions: 



* E-mail: izalamea@phys.columbia.edu 

t Also at Astro-Space Center of Lebedev Physical Institute, Prof- 
sojuznaja 84/32, Moscow 117810, Russia 



(i) The flow is axially symmetric. We assume that the 
rotational axes of the accretion flow and the black hole are 
aligned (or anti-aligned). 

(ii) The flow is symmetric under reflection about the 
equatorial plane. The symmetric streamlines collide in the 
equatorial plane and form a ring-like caustic around the 
black hole. 

(iii) The flow falls freely (ballistically) from a large ra- 
dius r ^ Tg until it hits the caustic; its pressure is negligi- 
ble everywhere except in the mini-disc. This assumption is 
valid if the flow is cooled efficiently (by radiation in X-ray 
binaries or by neutrinos in collapsars). The heat released in 
the shocks that accompany the disc-infall interaction is as- 
sumed to be quickly radiated away, so that the shocks stay 
near the disc plane, forming a "sandwich". The validity of 
this assumption is discussed in BIOl for the case of X-ray 
binaries and in Beloborodov (2008) for collapsars. 

(iv) The flow is quasi-steady: its accretion rate and an- 
gular momentum remain constant on the timescale of ac- 
cretion through the mini-disc (which is comparable to the 
free-fall time from r ~ lOr^). 

The paper is organized as follows. Section 2 describes 
the parabolic ballistic infall in Kerr metric and its ring caus- 
tic in the equatorial plane. In section 3 we write down the 
equations that govern the gas motion in the caustic (the 
mini-disc) and solve the equations numerically. In section 4 
we calculate the total luminosity of the disc observed at in- 
finity, taking into account the light capture into the black 
hole. 
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2 SUPERSONIC INFALL WITH ANGULAR 
MOMENTUM 

The spacetime of a black hole of mass M and angular 
momentum J is described by the Kerr metric. In Boyer- 
Lindquist coordinates {ct^r^O^cj)) the metric is given by 
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The spin parameter of the black hole a — J /Mc has dimen- 
sion of cm and must be in the interval \a\ ^ = GMjc? . 
Throughout the paper we shall also use the dimensionless 
parameter = a(? jGM^ |a^| ^ 1. 

We assume that the gas infall forms at a large radius 
r ^ rg, where it is efficiently cooled and, like dust, begins 
to fall freely towards the black hole. A streamline of this 
ballistic infall is determined by three integrals of motion: 
specific angular momentum /, its projection Iz on the spin 
axis of the black hole, and specific orbital energy E ^ (the 
infall is nearly parabolic). The four- velocity of a parabolic 
free-fall in Kerr metric is given by (e.g. Misner et al. 1973), 
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where h = /sin^oo, ^oo is the asymptotic polar angle of 
streamline at large r, and 
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The ballistic accretion flow is completely specified by the dis- 
tribution of its density and angular momentum on a sphere 
of a large radius r ^ rg. As the gas approaches the black 
hole and develops a significant rotational velocity, the infall 
is deflected from pure radial motion, and its streamlines in- 
tersect in the equatorial plane. The radius of this collision 
is determined by the angular momentum of the colliding 
symmetric streamlines. 



2.1 Collision Radius 

The streamline coming from an asymptotic direction 
(^oo,0oo) with angular momentum /(^oo,0oo) reaches the 
equatorial plane and collides with the symmetric streamline 
at the radius defined by 
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Figure 1. Fractional difference between the exact and its 
approximation (eq. Illj) as a function of I. The curves start at 
^0 = I* that corresponds to at the horizon of the black hole. 
All curves are plotted for maximally rotating black holes (a^ = 1), 
when eq. is least accurate. 
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We solved this equation for numerically (it involves an 
elliptic integral on the left-hand side). A good analytical 
approximation to is given by 
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2ac/sin^oo. The accuracy of this 
approximation is shown in Figure [T] It is better than 3 per 
cent when is outside the horizon, for any set of parameters 
relevant for the mini-disc formation. 

The collision radius is maximum for the streamlines 
near the equatorial plane (^oo 7r/2). This defines the outer 
radius of the disc, 



(12) 



2.2 Non-intersection of streamlines above the 
equatorial plane 

The ballistic flow defines a mapping from the sphere 
(^oo,0oo) to a sphere of radius r (cf. BIOl): (^oo,0oo) 
{0{r)^(t){r)). The streamlines of the flow do not intersect be- 
fore reaching the equatorial plane if the Jacobian of this 
mapping remains positive for all r > r^, i.e. 



/ de{r) d(i){r) \ 

oo 



J = Det 



V 



de{r) d<j){r) 



> 0. 



(13) 



This condition is equivalent to dO{r) / dOoo > 0. Using the 
relation between and r. 
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one finds that J' > if 
dcos 0{r) 



d cos Oo 
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This condition must be checked for a given distribution of 
angular momentum 1(0 oo) - In the numerical examples below 
we assume a distribution of the form 



1(00 
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(rigid-body rotation at infinity). We find that the stream- 
lines do not intersect before reaching the ring caustic in the 
equatorial plane. 



3 DISC DYNAMICS 

We consider here only accretion flows that are asymptot- 
ically spherical at r ^ r^, i.e. we assume that the accre- 
tion rate at infinity is spherically symmetric, dM/dQoo — 
const — Mtot/47r. Streamlines that start at ^oo with angular 
momentum l{0 oo) reach the equatorial caustic at a radius 
^^★(^oo). The accretion rate through the disc at a radius r, 
M(r), equals the net accretion rate along the streamlines 
that enter the caustic outside r, 



M(r) — Mtot cos^ 



(17) 



where Ooo{r) is the asymptotic polar angle of streamlines 
that collide at radius r. We find it by inverting the function 
r^(^oo). 

Matter inside the disc moves horizontally with four- 
velocity — {u* ^u^ ^ O^u^) and density p. The disc is steady 
and axially symmetric, so and p depend on r only. Equa- 
tions for u^(r) and p{r) are derived using the conservation 
laws for baryon number, energy and momentum. These laws 
are expressed by the following general equations (e.g. Lan- 
dau & Lifshitz 1980), 
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where T^^ = pc^u^v!' is the stress-energy tensor; we assume 
that it is dust-like everywhere, i.e. neglect the internal en- 
ergy density, pressure, and magnetic fields compared with 
pc^ . Above the disc, the stress-energy tensor is that of the 
ballistic infall. The infall four- velocity and density just above 
the disc shall be denoted by and p^ to distinguish them 
from the similar quantities inside the disc, and p. Then 
the conservation laws feas. [T8land [T9)) give. 
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Here M = 27v r h pu^ , h is the thickness of the disc, and 
T^^ = pc^u^u^ is the stress-energy tensor of the infall. Using 
equation (|20|) one finds from equation (|2T]) 
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For i — {t, 0} (conservation of energy and angular momen- 
tum) the right-hand side of this equation vanishes. In par- 
ticular, for i = this equation gives 
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For l{0oo) given by equation ([16} we have Ucj) = Iz = l^/lo 
and, using equation ()17|) . we find 
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For i = r, equation (|22]) expresses conservation of ra- 
dial momentum. Its right-hand side does not vanish, and 
one needs to evaluate ^^u^uK Substituting — {ucp — 
9tct)U^) / 94)4), using UiU^ = — and collecting terms one can 
write 
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Finally, can be expressed in terms of and from 
UiU^ — — c^. 
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In equation (|28p all quantities are known functions of radius 
except . We solve numerically this differential equation 
for u^(r). Example solutions are shown in Fig. [2l where we 
plot the radial and azimuthal velocities measured by ZAMO 
(zero- angular- momentum observer at fixed r; Appendix A 
gives the transformation of 4- vectors to the ZAMO frame). 

If the angular- momentum parameter of the flow, /o, ex- 
ceeds a critical value /cr, accretion in the disc is stopped 
by the centrifugal barrier {u^ changes sign). For flows with 
lo < lev the radial velocity remains everywhere negative. As 
/o approaches /cr the trajectory of disc accretion makes more 
turns around the black hole (see Fig. [3]) and at /o = Uv it 
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Figure 2. Upper panel: solutions for the radial velocity measured 
by ZAMO (in units of c) of the parabolic infall above the disc 
and the matter inside the disc. Two cases are shown: = 0.6, 
^0 = 2.1rgC and = 0.95, Iq = 1.56r^c. Lower panel: the corre- 
sponding azimuthal velocities of the infall and the disc, measured 
by ZAMO. 



makes infinite number of turns. Similar behavior was found 
for mini-discs around Schwarzschild black holes; in this case 
lev = 2.62r^c (BIOl). For spinning black holes, lev depends 
on a^. We have evaluated numerically /cr(a*) and the corre- 
sponding maximum size of the mini-disc rcz(a^). The results 
are shown in Figure 2] 

Figure U also shows the minimum value of /o, denoted 
by that is required to form a disc outside the black- hole 
horizon, 




Figure 3. Trajectory of disc matter in the critical case Iq = 
2.1ArgC for a black hole with = 0.6. The insert shows the 
radial and azimuthal velocities of the disc measured by ZAMO. 
X and y are the coordinates in the equatorial plane defined by 
X = r cos (f)^ y = r sin cf). 




i + \/r" 
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Accretion proceeds in the inviscid mini-disc regime when 
h < Iq < lev This range shrinks with increasing > 0, and 
expands if < (which means that the black hole and the 
accretion flow rotate in the opposite directions). 

Matter in the polar region of a quasi-spherical accre- 
tion flow falls directly into the black hole before crossing 
the equatorial plane. The mini-disc is formed in the other, 
equatorial part fo the flow (cf. Fig. 1 in BIOl). The boundary 
between these two accretion zones is determined by the con- 
dition r^(/) = Th- Equation (|17|) gives the fraction of Mtot 
that accretes through the disc. 
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which depends on Iq and a; it is maximum when Iq 
The maximum fraction is shown in Figure [5] 



Figure 4. Bottom panel: the range of angular momenta li,(ai,) < 
lo < ^cr(a^) that lead to the inviscid mini-disc regime (shaded re- 
gion). Accretion with < is quasi-spherical all the way into 
the black hole; it does not form a caustic outside the horizon. 
Accretion with Iq > ^cr must proceed through a viscous, centrifu- 
gally supported disc. Top panel: maximum radius of the mini-disc 
(/o = Icr) as a function of the black-hole spin parameter (solid 
curve). The radius of the black hole r/i(a^) (eq. I30|) is shown by 
the dotted curve. 



4 DISC LUMINOSITY 

The mini-disc is a radiative caustic that converts kinetic 
energy of the acretion flow to radiation. Let L(r) be the 
total luminosity produced by the disc outside radius r. The 
law of energy conservation gives an explicit expression for 
the luminosity. 



L(r) - [c^ + cut{r)] M{r), 



(32) 
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Figure 5. Maximum mass fraction accreted through the inviscid 
mini-disc. This maximum corresponds to Iq = lev and depends on 
the black-hole spin (see the text). 
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Figure 6. Radial distribution of escaping luminosity for discs 
with lo/rgC = 0.8, 1.1, and 1.5 around a black hole with spin 
parameter = 0.9. The luminosity is normalized by the rate of 
rest-mass accretion through the disc, M(r^)c^. 



where is the initial orbital energy of the accretion flow at 
infinity and —cut{r) is the orbital energy of the disc material 
at radius r. It is found using the solution for u'^{r) from 
section 3. The radial distribution of luminosity is given by 
dL/dr, and the total produced luminosity is L = L(r^). This 
luminosity would be received at infinity if no radiation were 
captured by the black hole. 

When the capture effect is taken into account, the lu- 
minosity escaping to infinity may be written as 



nr 
J ru 



dL 
dr 



/esc(r) dr, 



(33) 



where /esc(r) is the escaping fraction of radiation emitted 
at radius r; the fraction 1 — /esc is absorbed by the black 
hole. Assuming that the emission is approximately isotropic 
in the rest frame of the disc, we derive in Appendix A 
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where /3 is the disc velocity (in units of c) measured by 
ZAMO (zero-angular- momentum observer at fixed r), 7 = 
(1 - /3^)~^/^ and gtt = gtt - gtcj^/gcpcp- The integral is taken 
over the escape cone 5'esc(r) — all photon directions Q, (in 
the ZAMO frame) that lead to escape. The calculation of 
these cones is described in Appendix A. 

The resulting Lesc(^) is found numerically. Figure [6] 
shows the radial distribution dLesc/dr for mini-discs with 
three different lo around a black hole with spin parameter 
= 0.9. Most of the luminosity is produced in the region 
Tg < r < lAvg where the infall velocity relative to the disc 
is large and hence a large energy is released in the disc-infall 
interaction. 

We define the radiative efficiency of the disc as the ratio 
of the total luminosity radiated to infinity, Lesc = ^esc(^/i), 
to the rest-mass flux through the disc, 

Lesc rnr\ 

r]esc = ^— — — . (35) 
M(rh)c^ 

If the light capture into the black hole is ignored, i.e. Lesc is 
replaced by L, the efficiency is given by 



1 + 



ut(rh) 



(36) 



We evaluated numerically the dependence of 77 on the two 
parameters of the mini-disc Iq and (Fig. [7|). This depen- 
dence may be better understood if we express ut (vh) in terms 
of u^(rh) and U(f)(rh) from UiU^ — —(?\ then equation (|36p 
yields 



1- (1 + ^/13^)^^^ 



(37) 



If = this equation simplifies to 77 = 1 — \u^{rh)\/c 
and gives a monotonic dependence of 77 on Iq: flows with 
larger angular momenta have smaller \u'^{rh)\ (their radial 
motion is centrifugally decelerated) and higher 77. In the case 
of a rotating black hole, equation (|37|) has additional terms 
which lead to a complicated dependence of r] on lo and a^. 
For instance, when = 0.9 the dependence of 77 on lo is not 
monotonic (Fig.O. 

The efficiency 77 is generally increasing with increasing 
spin of the black hole. It is especially high for retrograde 
discs, which are described by the solutions with < 0. In 
this case, the second term on the right-hand side of equa- 
tion (|37|) is positive and can substantially increase 77. A re- 
markable feature of retrograde discs is that they can extract 
energy from the black hole via a hydrodynamic analog of 
Penrose process. This occurs if is close to —1; in such 
discs ut{rh) < and 7y > 1. The maximum 77 = 1.1 is reached 
when — —1 and lo — /cr(— 1) — 3.15r^c. 

Next, we evaluated numerically the efficiency 77esc that 
takes into account the capture of the produced radiation into 
the black hole (Fig. [8|). This effect greatly reduces the ob- 
served luminosity. The reduction is especially significant for 
retrograde discs because their radiation is Doppler-beamed 
in the direction opposite to the black hole rotation, so most 
of their radiation misses the escape cones shown in Fig- 
ure [AT] The resulting efficiency T^esc is highest for prograde 
discs around maximally rotating black holes (a^ = 1). The 
maximum 77esc is close to 10 per cent. 
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Figure 7. Radiative efficiency as a function of /q, ignoring the 
light capture into the black hole. Seven curves are plotted for 
black holes with different spin parameters a^. Solid curves are 
used for ^ and broken curves for < 0. 



to ~ YIGM I (?\ then it must switch to the standard viscous 
regime. 

We have calculated the radiative efficiency of mini-discs 
taking into account the effect of photon (or neutrino) cap- 
ture into the black hole (Fig. [5]). The efficiency is maximum 
when the mini-disc has its maximum size, near the transi- 
tion to the viscous-disc regime. For prograde discs around 
rapidly rotating black holes, the efficiency approaches 0.1, 
which is ^ 3 times higher than for non-rotating black holes. 

We have also studied the case of retrograde discs, where 
the disc and the black hole are counter-rotating (a^ < 0). 
Such discs can form in wind-fed X-ray binaries where gas 
accretes with alternating angular momentum due to fluc- 
tuations and 'flip-flop' instability (e.g. Shapiro & Lightman 
1976; Blondin & Pope 2009 and refs. therein). We find that a 
hydrodynamic analog of Penrose process works in such discs 
if is close to —1, i.e. the black hole is close to maximum 
rotation. Then the black hole accretes matter with nega- 
tive orbital energy, which means that energy is extracted 
from the black hole. However, if the extracted energy is ra- 
diated quasi-isotropically in the disc rest-frame, most of the 
produced radiation ends up inside the black hole, and the 
escaping luminosity is suppressed. 



REFERENCES 




u ^ — ' ^-r-^ — ^ — ^ — ^ — I — ^ — ^ — ^ — ^ — I — ^ — ^ — ^ — ^ — I — ^ — ^ — ^ — ^ — I — ^ — ^ — I 

0.5 1 1.5 2 2.5 3 

Figure 8. Same as Fig.[7]but taking into account the suppression 
of radiative efficiency due to light capture into the black hole. 



5 CONCLUSIONS 
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APPENDIX A: CALCULATION OF THE 
ESCAPING LUMINOSITY 

Al Angular distribution of emission in the local 
ZAMO frame 

Consider an infinitesimal element of the disc, a ring of radius 
r and thickness 5r. It produces luminosity 

SL=^6r= 4- lie" + cut)M] Sr. (Al) 
dr dr ^ 

The ring emits photons in all directions, which will be pa- 
rameterized by unit 3D vector Cl in the local frame of ZAMO 
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(zero- angular- momentum observer at fixed r, see e.g. Misner 
et al. 1973). Then 6L may be written as 



6L-. 



[ mdn= [ 



dt E djSL) 
dtE dQ 



(A2) 



Here 6L is the luminosity of the ring measured by ZAMO 
and t is the proper time of ZAMO. E / E is the ratio of photon 
energy measured at infinity, E = —cvt, and measured by 
ZAMO, E = —cvt = cv^] this ratio depends on the emission 
direction ft. It is found from the Lorentz transformation 
of the photon 4-velocity from the coordinate basis to the 
(orthonormal) ZAMO basis. 



V-9tt 







V — y^QrrV , 

where gtt - gtt 



/goov , 

^. Then one finds. 



E 
E 



Vt_ 

Vt 



-gtt 



jt^^i 



(A3) 



(A4) 



(A5) 



The ratio di/dt appearing in equation ()A2|) equals [—guY^'^. 

The angular distribution of luminosity measured by 
ZAMO, d(SL)/dQ, is related to the angular distribution of 
luminosity in the rest frame of the accreting gas of the disc, 
d(6Lc)/dQcj by the Doppler transformation (see e.g. Rybicki 
& Lightman 1979), 



d{6L) 
~dQ~ 



d{6Lc 



74(l-/3•^^)3 dQc 



(A6) 



where p is the disc velocity (in units of c) measured by 
ZAMO. We assume that emission is approximately isotropic 
in the gas frame, i.e. d{6Lc)/dQc = 6Lc/4:7v. Substitution of 
equations ()A5|) and (|A6p to equation (|A2p gives, 



6Lc , ~ \i/2 



{\J-~gtt - {gt4>l^g^)^^\ 
74(1-/3-1^)3 



dVi. (A7) 



The expression for escaping luminosity (5Lesc is similar ex- 
cept that the integral is taken over the escape cone ^Sesc 
rather than 47r. Therefore, the escaping fraction is given by 



SLe 



6L 



W-gtt - {gtcf>/^/g^)^'^] ^ 



W-gtt - {gtct>/^/g^W] 
74(1-/3- A)3 



Evaluating the integral over 47r we obtain equation (|34p . 



Using the relations between and v"^ feas. lASl and IA4|) . we 
express Cl in terms of t'*. 



sm a sm — 



Qrr _ Y^^^^^+~^M^^+^^ 

gtt + r(r — Vg) ' 

_gee_v^ _ 
gtt V* 



+ a'^r{r + rg) v 

q2 _|_ _ j.^^ yt ' 



(All) 



(A12) 



COS a = 



where all metric coefficients have been evaluated at the equa- 
torial plane at the emission radius r. 

The photon 4-velocity is expressible in terms of four 
integrals of motion in Kerr metric: energy E, total angular 
momentum L, its projection L^, and Carter integral F. Using 
these expressions (see e.g. Chandrasekar 1982) and choosing 
the affine parameter along the photon worldline so that E — 
c^, we find /v^ and /v^ as functions of r, F and L^, and 
obtain 

+ a^r(r + Vg) 
[r^ + a'^r{r + Vg) — TgaLzv]^ 
X [r^ + (F + a") + rvg {a -T - 2aL,) 

W{T + lI)], (A13) 



sin^ a sin^ if — — {V + L^) (r^ + — r^r) 

[r^ + a^r{r + rg)\ 

[r^ + a^r{r + r^^) — TgaLzv] 

Equations ()A13|) and ()A14|) can be solved for Lz and F for 
given r, a and Lp. 

The fate of a photon is determined by the equation of 
radial motion for null geodesies (e.g. Chandrasekar 1982), 

[r^ + cos^ {^^Y — (^^ + — r^r) {V — a^E^^ 

+ (r^ + a^)^ - 2rgaELzr + a^L^, (A15) 

where E = with our choice of affine parameter along pho- 
ton worldline. The possible turning points rtum are found 
from the condition = 0. Let ro be the initial radial po- 
sition of the photon. Two cases are possible: (i) v{ro) > 0, 
the photon escapes if there are no rtum > tq. (ii) v(ro) < 0, 
the photon escapes if there is at least one turning point such 
that rh < rturn < ro and no rtum > ro- 

For any given direction ft, we find Lz and F from equa- 
tions ()A13P and (|A14p . then determine the roots rtum of 
equation v^{r) — and check the escape conditions. All es- 
caping directions form the "cone" Sesc on the sky of ZAMO, 
which is found numerically. Figure lAll shows the escape 
cones for five emission radii ro and four values of a^. 



A2 Escape cones 

The photon direction in the ZAMO frame, Cl, can be de- 
scribed by two angles a and cp, 

—r 

or = cos a, (A8) 

yt 

-6 

Cl^ — — — sinasin(/?, (A9) 

yt 

^ = —- — sinacos(/?. (AlO) 

yt 
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Figure Al. Escape cones Sesc on the ZAMO sky. The photon 
direction Cl is specified by two angles a and feqs. lASIAlOj) . The 
figure uses polar coordinates {p,(f)) with p = a/ir to represent 
all possible photon directions. The origin of the diagram p = 
corresponds to the radial direction away from the black hole 
(such photons always escape) and the unit circle p = 1 (thick 
black curve) corresponds to the radial direction into the black 
hole (such photons are captured). The colour curves show the 
boundary of the escape cone for five emission radii; the emission 
radius is indicated next to the curves, in units of the horizon 
radius r^. The figure presents four such diagrams calculated for 
black holes with spin parameter = 0, 0.6, 0.9, and 0.998. 



